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Abstract 

Cosmic ray signals from dipole-interacting dark matter annihilation are considered in the positron, 
antiproton and photon channels. The predicted signals in the positron channel could nicely account 
for the excess of positron fraction from Fermi LAT, PAMELA, HEAT and AMS-01 experiments 
for dark matter mass larger than 100 GeV with a boost (enhancement) factor of 30 — 80. No excess 
of antiproton over proton ratio at the experiments also gives a severe restriction for this scenario. 
The predicted signals as monoenergetic gamma-ray lines (monochromatic photons) for the region 
close to the Galactic center are investigated, and we show that they are in the current observational 
experimental constraint. To the end, the gamma-ray excess of recent tentative analyses based on 
Fermi LAT data and the potential probe of the monochromatic lines at a planned experiment, 
AMS-02, are also considered. 
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I. INTRODUCTION 



The existence of the dark matter (DM), as the invisible matter interacting by the force of 
gravity, has been widely accepted by cosmological observations from the experiments: Cos- 
mic Microwave Background (CMB) Galactic Rotation Curves [2], Gravitational Lensing 
[jj and Massive Compact Halo Objects and etc. About 83% of the matter (around 23% 
of the total energy density) in the universe is believed to be composed of the DM to account 
for the observations. However, the nature of DM is still completely unknown despite decades 
of detection efforts. Many possible explanations have been proposed. One of the alternative 
explanations from the point of view of particle physics is that the DM is composed of massive 
particles and its interaction with ordinary matter is very weak. 

Recently several DM models 1 (inelastic DM [7, 8], asymmetric DM {9], form factor DM 



10| ) by magnetic dipole interaction have been considered. The DM in these models have 



a few states and have no direct interaction with the photon. The candidate particles could 
thus be stable and make up the invisible matter in our universe. On the other hand, the 
direct interaction with photon through magnetic dipole coupling has gotten some more at- 
tention 
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17] due to its plausibility. Most of works concentrated on direct detections of the 



DM. The main motivation of magnetic dipole interacting DM scenario is that the magnetic 
dipole coupling can be sizable compared to other electromagnetic couplings, because the 
magnetic dipole conserves the discrete symmetries like parity (P), time reversal (T), and 
charge conjugation (C) or its combinations. 

In this work, we consider cosmic ray signatures (indirect detections) of the direct dipole- 
interacting DM with the shifted photon (hypercharge gauge boson B^). The cosmic ray 
signatures in the positron, anti-proton and photon channels are considered for the DM mass 
near the electroweak scale (10 — 1000 GeV), essentially around 100 GeV. The dimension-5 
operator which induces the dipole interaction is ifja^ipB^ , and it may be expressed about 
photon and Z boson in the standard model context since hypercharge gauge boson is a 
linear combination of photon and Z boson with Weinberg angle 9w- The relevant effective 



1 Actually these models have been built to explain the annual modulation signal from DAMA/Nal Q 
and DAMA/LIBRA Q experiments with null results from other experiments (there is no experimental 
evidence corroborating this signal yet). The scenario in Ref. [7[ is especially interesting, because the signal 
appears electromagnetic energy deposit, not nuclear energy deposit, through the single photon emission 
by the decay of the excited state. 



2 



Lagrangian is given by 

C cS = ^a^(F^ - taaOwZ'"), (1) 

where F^ v is the field strength for photon, Z^ u for Z boson and \x is the magnetic dipole 
moment. The DM annihilations therefore proceed to the standard model particles via 7, Z 
exchanges. The annihilation processes were studied in Ref. [ll| ii 



in detail, and here we take 



an advantage of the results (annihilation rates or branching ratios). 

II. COSMIC RAY SIGNATURES 

The DM may annihilate at some point in Galaxy and produce the standard model parti- 
cles. These produced particles then propagates in the interstellar medium. The anti-matters 
and photons have been considered to be the subject of indirect DM searches, because the 
anti-matters are rarely produced in astrophysical process and gamma rays can transport 
freely without energy loss or transmutation of the direction. They may thus provide impor- 
tant signatures of the DM in the Galaxy. Observations of such signals can reveal information 
on the microscopic nature of the DM. 

The emissivity /energy (production rate or source for the signals) 2 at location x from the 
Galactic center is obtained from the convolution over the various annihilation channels / of 
the annihilation rate {ov) ^ with the differential yield (single particle spectra) 3 (jdN^/dTj 
for the final state particles a, 

where M is the DM mass, B is an overall boost (enhancement) factor and p(x) is the 
DM mass density at the location x. The DM mass density about the location from Galactic 
center (DM halo profile) is not known, especially nearby center (< 100 pc). The theoreti cally 



motivated ones are Navarro- Frenk White (NFW) 23j, Moore 24|, cored Isothermal 25], 



2 If the DM is produced with a primordial asymmetry like baryons, there would be almost no signals from 
DM-antiDM annihilations due to lack of anti-darkmatters (antiDMs). This source is for equal populations 
of the DM and antiDM, and hence our predicted signals will be upper limits of the predictions. Recently 
a mechanism of the DM-antiDM oscillations is suggested to re-equilibrate the populations at late times 
00. 



3 We use PYTHIA [20|, as implemented in DarkSUSY [2l| or MicrOMEGAs [2g| program, to generate the 
differential yields (injected particle spectra). 
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and recently Einasto profiles [2g,|27(|. The kinetic energy is often approximated to the total 
energy E, in case that the particles are energetic. Notice that the factor | is different from 
the one of self annihilating DM. 

The charged particles which are generated from the DM annihilation are predicted to 
come from the halo near the Sun, not too far from the Sun at least, because they may 
lose the energy while propagating through the Galactic halo. They are deflected by the 
Galactic magnetic field, and this property has been described by space diffusion J^. The 
charged particles also suffer energy losses from synchrotron radiation and inverse Compton 
scattering. Solar modulation can induce a certain amount of energy losses. Their energy 
spectrum at the Earth, therefore, differs from the one produced at the source. The equation 
that describes the evolution of the energy distribution for the charged particles may be 
expressed as 

= Q a (x,T)-2hT ann 5(z)^—j , (3) 

where dn/dT is the number density of particles per unit volume and energy. The second term 
accounts for the space diffusion with the energy dependent diffusion constant K (T). The 
energy loss due to synchrotron radiation and inverse Compton scattering on CMB photons 
and on Galactic starlight is described in the third term with the rate of energy loss b (T), and 
the fourth term deals convective wind. The last term accounts for the annihilation of the 
produced matters in the interstellar medium, H and He atoms, with annihilation rate T ann 
in the disk of thickness 2h ~ 0.2 kpc, and hence it is provided as a negative source term. 
The relevant coefficients were parameterized, and the established (transport) parameters 
have been estimated from the analysis of observed isotope ratios in cosmic rays, primarily 



the boron to carbon (B/C) ratio 29]. Three propagation models have been featured with 



the established parameters, and these propagation models correspond to minimal (MIN), 



medium (MED) and maximum (MAX) anti-proton fluxes 29. l30j|. 



The number density dn/dT is obtained by solving Eq. (3) with the steady state condition 



~§t (Jr) = an d boundary conditions in a two- zone model 3l|, where the region of diffusion 
of cosmic rays is represented by a thick disk of thickness 2L ~ 5 — 20 kpc and radius R ~ 20 
kpc, and the thin Galactic disk lies in the middle and has a thickness 2h, radius R. The 
boundary conditions are such that the number density vanishes at z — ±L and at r = R. 



A. The Positron Channels 



The propagation of positrons and electrons is dominated by space diffusion and energy 
losses 

(4) 

This diffusion equation may be solved by the Green function formalism or the Bessel- 
transform method, and the solution results in the following form, 

■■M 



dn \ 1 f M 

dE) e+ = b{E)j E dEsQe+ (X0 ' Es) Ge+ (E ' Es 



(5) 



where x Q is the location of the Sun from the Galactic center. The function G e +(E, Eg) must 
fully encode the Galactic astrophysics from the input e nerg y Eg to energy E(< Eg), and 



the full expression of this function can be found in Ref. 30]. 
The positron flux is then given by 



.v P + 



dn \ 



= v+(<n,), ,PoV f M (dNl\ 



where p ~ 0.4 GeV/cm 3 32[ is the DM density in the solar vicinity, and v e + is the velocity 
of the positron. 

The enhancement may come from subhalo structure (dark clumps), and the non- 
perturbative Sommerfeld effect because one of force carriers is photon in this scenario. The 
enhancement by the Sommerfeld effect can be calculated from the original form in this case. 
We can split the dipole operator in energy and momentum dependent parts by the familiar 
Gordon decomposition, 

/JAJ(p')a^q u u(p) = i/jJU(p)(2M-f fl + (p' - pY)u{p). (7) 

In the energy dependent part, we have the same type of coupling with the DM mass depen- 
dence as electric charge. The original Sommerfeld enhancement factor is 

where v is the DM velocity (~ 10~ 3 ). In the original form, a M is the electric fine structure 
constant, but a M = lGnfi 2 M 2 in this case. We have the enhancement S ~ 16 for the 
DM mass of 100 GeV and magnetic dipole /i ~ 0.1 TeV _1 (~ 10~ 17 e-cm). The strength of 



magnetic dipole is, according to Ref. 



Ill ], almost constant for the DM mass larger than 100 




FIG. 1: Positron fraction from the annihilation of the fermionic dark matter particle. The boost 
factors have been chosen to provide qualitatively good fits to the data with dark matter mass 100 



GeV. Shown are the background in solid-line and the experimental datasets of the PAMELA 
HEAT 36], AMS-01 37], Fermi LAT [38], and CAPRICES94 3]- 



ad], 



GeV, and thereby the enhancement increases in the square of the DM mass. Thanks to this 
reason, the fluxes, Eq.(6), must have the similar magnitude about the DM masses larger 
than 100 GeV because the fluxes are not scaled by the DM mass in our scenario. We notice 
that the fluxes are scaled by the DM mass in most of models, and large boost factors have 
been required for larger DM mass because the fluxes decrease in the square of the DM mass. 

The positrons are affected by solar wind and lose energy while transporting in the so- 
lar system. This effect leads to a shift in the energy distribution between the interstellar 
spectrum (IS) and the spectrum at the top of the atmosphere (TOA). This modulation is 
considered for the predicted fluxes by the relation in force field approximation [33j, 



E 2 



m; 



as 



{E+ \Z\e$) -m 2 + 



(E+\Z\e$), 



(9) 
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where \Z\ is the magnitude of electric charge (1 in this case), e is the electric constant and 
$ is the Fisk potential, namely solar modulation parameter, which varies between 500 MeV 
and 1.2 GeV over the eleven- year solar cycle. Since experiments are usually undertaken 
near solar minimum activity, we choose $ = 600 MeV (the Fisk potential for the PAMELA 
experiment) for our numerical analysis. 

We show in Fig. 1 the predicted positron fractions for the boost factors, B = 30, 50, 80, 



with the computed background from Ref. 



34| and several experimental data sets. The 



predicted fluxes have almost no difference in the halo profiles or the diffusion models. The 
predicted fractions exhibit a rather sharp distribution at E e + ~ M, since our candidate can 
directly annihilate into electron and positron pair. The PAMELA 35j have shown a steep 
rise in the 10 — 100 GeV range in their measurements and confirmed the results of HEAT 



361 ] and AMS-01 [37J experiments. Recently, the steep rise has been extended to 200 GeV 



with three more data points over 100 GeV at the Fermi LAT 38j. The predicted signals 
with the boost factor 30 — 80 nicely fit measurements of the PAMELA for the DM mass of 
around 100 GeV. An enhancement of about a factor of 16 comes from the Sommerfeld effect, 
and the rest, an enhancement factor of 2 — 5, is expected to come from subhalo structure 



(dark clumps). The typical subhalo boost factor has been predicted to be ten 40j, |4l|] or 
twenty [2] at most. Our subhalo boost factors are in the reasonable range. We also have 
the difference between predictions and experimental measurements or background at low 
energies (< 10 GeV). It has been noticed that solar modulation effect we consider has no 
charge-sign dependence and has to be modified. This must be a future study. 



B. The Anti-proton Channels 

The propagation of anti-protons is dominated by diffusion and the effect of the Galactic 
wind. The energy spectrum of anti-protons is obtained by solving the following diffusion 
equation, 

|(^ W V O (gJ-V.( K mvg)J=g.(x,7 1 -2 H -„i W (g / (10) 

An important difference with the positron case is that energy loss of anti-protons is negligible, 
because anti-protons are more massive and hence it is absent in the diffusion equation (10). 
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FIG. 2: Antiproton over proton ratio as a function of kinetic energy. Shown are the background in 

n 

CAPRICE98 [48] and 



solid-line and the experimental data sets of the PAMELA 

IMAX y. 



46], BESS 



1UW I J 



The anti-proton flux is then given by 



4tt \dT 



B 



v ¥ (av) 



16n 



(dNf\ / 



PQ 

M 



GAT) 



(11) 



p \ / p 

where the function Gp (T) encodes all the astrophysics, and the full expression can be found 



in Refs. 



43 



Fig. 2 shows the predicted ratios of anti-proton over proton with the computed back- 
ground [45| and several experimental measurements. As in the case of the positrons, the 
predictions have almost no difference in the halo profiles, but it is sensitive to the prop- 
agation models. The predictions in the MIN propagation model are selected, and they 
are within the unobservational range. However, this scenario is most likely ruled out for 
other propagation models 4 , MED and MAX, because the predicted fluxes in MED or MAX 



To be strict, the MED propagation model can also be viable because of the uncertainty of the transport 



S 



propagation models are about ten or one hundred times larger than the ones in the MIN 
propagation models. 



C. The Gamma Ray Channels 

The production of gamma-rays has been considered to be a very important channel to 
search for the DM signals, since they travel in straight lines and can travel greater distances 
without energy loss. For these reasons, they contain spectral and directional information 
that can be well measured. The gamma-ray flux from the DM annihilations at a given 
photon energy from a direction that forms an angle ip between the direction of the Galactic 
center and that of observation is accounted for by the line-of-sight (los) integration method, 



H/ fdNf\ f 2f , m , 



^ B ^A^)X p{rM))ds ' (12) 

where r = + s 2 — 2r s cosip is the Galactocentric distance. In terms of the galactic 
latitude b and longitude I, one has cos^ = cos b cos/. The coordinate s parameterizes the 
distance from the Sun along the line-of-sight. 
This form can be reduced to 

n , , n-1 B(av) f /100 GeV\ 2 /dJVA - , . 

8.31X10- (cm 2 - sr. s) . yF ^ T {-^ ) • ^ , (13) 

where Jac is a dimensionless line-of-sight integral averaged over the solid angle AQ and is 
defined by 



5 7 ~ 



with 

JM = /^p (15) 

J Pe r & 

In addition to the continuum emission, the direct DM annihilations produces 77 and 7Z 
final states 5 in this scenario, in which Feynman diagrams are shown in Fig. 3. Such processes 
would yield the very distinctive feature of monoenergetic gamma-ray lines (monochromatic 

parameters. The parameters are assigned to the best fit of cosmic ray B/C data [2^. However, the 
assigned values of transport parameters may differ by one order of magnitude. 
5 The production of the single 7 or Z final state is prohibited, because this process is impossible to conserve 
energy and momentum together. 
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FIG. 3: Feynman diagrams for annihilation of dark matter to one-photon or two-photons. The 
hatched circles indicate the dipole couplings. 



photons) with an energy E 1 = M or M{1 -m|/(2M) 2 ). The full annihilation rates for 

> 77 = /i 4 M 2 /2vr and (av) z 



the production of 77 are (<jv) = \i M 2 /In and (av) z = /z 4 M 2 /3| tan 2 6w/^ for jZ. The 



contributions for 77 and 7Z final states are suppressed by the magnetic dipole //, and the 
jZ final state has an additional suppression with Weinberg angle tan 2 6w- 

We restrict our analysis to the possible signals from the Galactic halo as in the case of 
the positrons for a complementarity, and hence the same boost factors 6 B = 30, 50, 80 are 
selected. We compare our predictions to experimental observations in two stringent cases. 
One is a process that may contribute to the extragalactic gamma-ray background (EGB). 
The other is a process to satisfy the experimental exclusion limit of the Fermi LAT [50J, in 
the case that astrophysical sources account for the EGB in entire energy range. 

Fig. 4 shows the predicted gamma-ray spectra as a function of photon energy in the region 
0° < £ < 360°, |6| > 10°. The spectra are superimposes of the continuum and monoenergetic 
gamma-rays at the DM mass of 100 GeV. The spectra have almost no difference in the halo 



profiles. The EGB from the Fermi LAT 50J is given by 



, E n-0.41 

£ 2 7 ~ 5.5 x 10" 7 f J (cm 2 • sr • s) GeV. (16) 

The same type of background from the analysis of the EGRET measurements is also de- 



scribed [5JJ with the one from the first analysis [52]. The predicted spectra are within the 
EGB, and we might have a signature if it can be disentangled from astrophysical ones. In 
addition, we check if our prediction can account for the EGRET anomaly which is not con- 



6 The enhancements can be different in each channel on the subhalo structure. Here we select the same 
boost factors in the assumption that the enhancements by the subhalo structure are not too different in 
each channel. 
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FIG. 4: The gamma-ray spectrum from Galatic halo for dark matter mass of 100 GeV and the 
region, 0° < t < 360°, \ b\ > 10°. Shown are the experimental data sets from the Fermi LAT and 
EGRET with their fitted spectral indices. 

firmed at the Fermi LAT. Our predictions are softer to explain the observation of the EGB, 
and the more enhancement would be needed. 

The predicted fluxes have to be within the uncertainty of Fermi LAT data, in the case 
that the EGB is accounted for by astrophysical sources in entire energy range. Fig. 5 shows 
the predicted fluxes with 90%, 95% and 99% experimental exclusion limits of the Fermi LAT, 
which are estimated from the data table in Ref. (so | . The predictions are under exclusion 
limits and are satisfied with the current unobservable experimental constraint. 

The Galactic center or the region close to that must be the most complex regions in the 
Galaxy due to many possible sources and the difficulty to model the diffuse emission. Hence, 
it may be very difficult to disentangle possible DM annihilation signals from the background 
fluxes. The monochromatic gamma-ray lines appeared from DM annihilations could provide 
smoking-gun signatures for these regions, because the line signals cannot be mistaken for 
astrophysical source. 

Fig. 6 shows the predicted annihilation rate for 77 and 7Z final states as a function of the 
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FIG. 5: The gamma-ray flux from Galatic halo for dark matter mass of 100 GeV and the target 
region 0° < i < 360°, \b\ > 10° with 90%, 95% and 99% experimental exclusion limits of the Fermi 
LAT. 



DM mass. The Sommerfeld effect is included, but not the subhalo structure. The strength 
of DM magnetic dipole is chosen to satisfy the relic density from Ref. U|. The curves 
given for the NFW, Einasto, and isothermal DM distributions are 95% CL. upper limits of 
Fermi LAT [53( for the region \b\ > 10° plus a 20° x 20° centered at the Galactic center. 
The predicted annihilation rates are approximately 4. Ox 10 -29 cm 3 /s at E 1 ~ 100 GeV for 
the 77 final state and 4.0 x 10~ 30 cm 3 /s for the 7Z final state at E 1 ~ 80 GeV. Otherwise, 
the estimated upper limits at the Fermi LAT are 10~ 27 — 10~ 26 cm 3 /s for both final states, 
depending on the DM mass and the halo profiles. The predictions are two or three orders 
smaller than the experimental upper limits. However, our predictions can be enhanced 
if dense DM clumps are considered in regions close to the Galactic center. Recently, the 



authors in Refs. 



54 



55J pointed out the gamma-ray excess, 1 — 3 x 10 cm /s, around 130 



GeV in the spectrum based on these measurements with 4.5 or 6a statistical significance. 
Our predictions with the subhalo boost factor of about 100 can account for the gamma-ray 
excess. 
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Dark Matter Mass (GeV) Dark Matter Mass (GeV) 



FIG. 6: Annihilation rates with Sommerfeld enhancement factor S for dark matter annihilation to 
77 or jZ. Shown are the 95% C.L. upper limits of the Fermi LAT for the region \b\ > 10° plus a 
20° x 20° centered at the Galactic center. 



One of the reasons we mostly get the low predictions for sharp peaks is due to the relatively 
poor energy resolution. The current energy resolution of the Fermi LAT is 11 — 13% [53J for 
the full width at half maximum. The energy resolution could be as better as 1.5 — 2% for 
a planned experiment, AMS-02 56|, in which the upper limits of annihilation rate will be 
10~ 30 — 10 -29 cm 3 /s. Our predicted signals are in the potential probe at the AMS-02. 



III. CONCLUSION 



We considered cosmic ray signals in the positron, anti-proton and photon channels of 
dipole-interacting DM annihilation. The predicted signals in the positron channel could 
nicely account for the excess of positron fraction from Fermi LAT, PAMELA, HEAT and 
AMS-01 experiments for the DM mass larger than 100 GeV with a boost factor of 30 — 80. 
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We estimated that an enhancement of about a factor of 16 would come from the Sommerfeld 
non-perturbative effect and the rest, an enhancement factor of 2 — 5, from subhalo structure 
(dark clumps). The predicted signals have almost no dependence on the DM mass because 
of the Sommerfeld effect. No excess of anti-proton over proton ratio at the experiments also 
gives a severe restriction for our scenario. This scenario must be viable for MIN Galactic 
propagation model, but likely ruled out for other propagation models, MED and MAX. 
The predicted signals from the Galactic halo in the region 0° < I < 360°, \b\ > 10°, and 
signals as the monoenergetic gamma-ray lines (monochromatic photons) for the region (\b\ > 
10° plus a 20° x 20° centered at the Galactic center) close to Galactic center were also 
considered. The predicted signals from the Galactic halo must be in the current unobservable 
experimental constraint, and the signals for the region near Galactic center as monoenergetic 

AT. The gamma-ray excess 



5i 



55] pointed out, could be 



lines are under experimental exclusion limits of the Fermi 
1 — 3 x 10~ 27 cm 3 /s around 130 GeV, which authors in Refs. 
accounted for in this scenario with the subhalo boost factor of about 100. Our predicted 
signals as monoenergetic lines for the region near Galactic center are also in the potential 
probe at a planned experiment, AMS-02, with the better experimental method. 
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